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We propose and study a realistic model for the decoherence of topological qubits, based on Ma- 
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fluctuating charge on a capacitively coupled gate, modeled by non-interacting electrons. In this 
context, we clarify the role of quantum fluctuations and thermal fluctuations and find that quantum 
fluctuations do not lead to decoherence, while thermal fluctuations do. We explicitly calculate decay 
times due to thermal noise and give conditions for the gap size in the topological superconductor 
and the gate temperature. Based on this result, we provide simple rules for gate geometries and 
materials optimized for reducing the negative effect of thermal charge fluctuations on the gate. 
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I. INTRODUCTION 

The central idea of topological quantum computing 
(TQC) 1-5 is to encode quantum information in a degener- 
ate ground state. In this context, the word "topological" 
means that one must use a system in which the ground 
state degeneracy cannot be lifted by any local and suffi- 
ciently weak perturbation^ Such a system is said to be 
topologically nontrivial. If it is prepared in one of its 
many possible ground states, it will stay in exactly this 
state as long as the system is only weakly perturbed. In 
other words, the information encoded in this way is free 
of decoherence. The manipulations of the quantum in- 
formation are performed via relatively large changes of 
material parameters that may be controlled externally, 
for instance via voltages applied to electric gates. The 
parameter variations needed in order to perform a non- 
trivial operation are so huge that they practically cannot 
happen accidentally. Moreover, the manipulation does 
not depend on the details of the parameter changes, but 
only on their topological properties. Thus, small imper- 
fections in the parameter changes are not supposed to 
lead to imperfections in the operation on the quantum 
information. All this is true in an ideal world in which 
the system stays in its degenerate ground state and never 
gets excited. 

One realization of such a topologically nontrivial sys- 
tem is the topological superconductor (TSC)^^— . As in 
any superconductor, a finite energy > A must be paid 
in order to excite the system above its ground state. 
However, the ground state itself can be degenerate in 
a TSC; in addition to the fcrmionic excitations with en- 
ergies > A, there are M "excitations" with zero energy, 
with M the number of pairs of spatially separated Ma- 
jorana bound states (MBS). Thus, the ground state is 
2 M -fold degenerate. For nontrivial rotations in this de- 
generate ground state space, the localized MBS must be 
moved (braided) around each other; such a braiding re- 
quires to change the local material parameters by large 



amounts, large enough that they cannot be due to un- 
controllable fluctuations. 

Recently, the concept of decoherence-free TQC, based 
on Majorana bound states in TSCs, has been challenged 
by the authors of Refs. [13 and[ll]. In these works it was 
emphasized that topological protection does not protect 
against all kinds of perturbations that might be present 
in a realistic experiment. In particular, the work of Bu- 
dich et aZ»ii investigates tunnel couplings of the Majo- 
rana bound states to nearby fermionic baths that are not 
gapped. Since the topological protection of Majorana 
based qubits is based on the conservation of fermionic 
parity (moving one single electron in a superconductor 
costs energy > A), but tunneling between ungapped 
fermionic baths and Majorana states allows single elec- 
trons to enter the TSC at zero energy, the quantum in- 
formation in such a system is lost on exponentially short 
timcscalcs. Thus, non-superconducting tunnel contacts 
close to the Majorana bound states must be avoided. 

Goldstein and diamond on the other hand, consider 
gapped fermionic baths which are coupled to the Majo- 
rana bound states by a fluctuating bosonic field, such as, 
for instance, phonons. In their work it becomes clear that 
a decay of the information (encoded by the local fcrmion 
parity) into a fcrmionic bath with an excitation gap A 
can occur only if the spectrum of the bosonic field is non- 
zero at frequencies u) ~ A. In other words, the energy 
that is needed to excite the TSC must be provided by 
the environment. One can also express this statement in 
terms of the adiabaticity requirement of TQC: if the ex- 
ternal parameters of the Hamiltonian describing the TSC 
(e.g., electric potentials) arc changed non-adiabatically, 
i.e., faster than H/A, then the TSC will not remain in its 
ground state, which is, however, an inevitable require- 
ment for TQC. 

Thus, it is clear that in the real world, where the con- 
stituents of an experimental setup are subject to all sorts 
of random fluctuations, the TSC will eventually get ex- 
cited and lose its quantum information. Thus, an explicit 
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study of the effect of different kinds of fluctuations on the 
TSC is inevitable for judging the viability of the concept 
of TQC. In this context, Ref. US raised the issue of deco- 
herence caused by quantum fluctuations of bosonic fields. 
In some instances, this work is interpreted to the effect 
that a bosonic quantum bath at zero temperature may 
cause decoherence of a topological qubit. It is one of the 
goals of the present work to clarify the role of quantum 
fluctuations, or, more explicitly, zero temperature fluctu- 
ations. 

Another much more obvious source of decoherence is 
the thermal fluctuation of a bosonic field. It is usually 
stated that the problems caused by those thermal fluctu- 
ations are exponentially small ~ exp(— A/T). However, 
the exponential dependence alone is not enough for judg- 
ing the feasibility of the concept of TQC - the pre-factor 
also matters. To be concrete, a decay time of a femtosec- 
ond times an exponentially large factor might still be ex- 
tremely short for the experimentally accessible range of 
temperatures. Thus, it is of utmost importance to actu- 
ally calculate the decay time of a topological qubit. Gold- 
stein and Chamon made a first step in this direction^ 
but they considered a generic model. Thus, no predic- 
tions could be made about the actual decay time but 
only about its functional form. In order to make explicit 
predictions about the feasibility, one must leave the level 
of generic models and focus on a particular experimental 
setup. This is the path we will go in this paper. 

The present work is based on an explicit experimental 
setup and a realistic model for one particular decoher- 
ence process, which we think is least avoidable in this 
setup (apart from quasiparticle poisoning, see Refs. [l2l — 
Il8| ). This explicitness finally allows us to plug in num- 
bers and to obtain bounds for various parameters, e.g., 
the temperature or the superconducting gap. The ex- 
perimental setup we choose has been proposed by Alicca 
and coworkers in Ref. Il9l . namely the movable Majorana 
bound states in networks of TSC wires. One of the ba- 
sic building blocks of this proposal is an array of electric 
gates by which individual sections in the wire network are 
tuned into or out of their topologically nontrivial phases. 
The Majorana bound states live at the phase boundaries 
so that a MBS can be moved by means of changing gate 
voltages. In Ref. [l9l these gates were assumed to be classi- 
cal objects giving rise to non-fluctuating local potentials 
that can be changed as slowly as needed. We extend this 
model by allowing charge fluctuations on the gates and 
study their effect on the TSC. 

The paper is organized as follows. In Sec. II we dis- 
cuss the setup and basic concepts. In Sec. Ill we derive 
an extended model for a TSC, including charge fluctu- 
ations on a nearby gate. The resulting decay time of a 
Majorana-based qubit is calculated in Sec. IV. In Sec. V 
we discuss our results and finally conclude in Sec VI. 



II. PRELIMINARIES, SETUP AND 
ASSUMPTIONS 

In order to give a self-contained and pedagogical de- 
scription of the decoherence process and to further mo- 
tivate our analysis, we start with some general consider- 
ations. 



A. Local fermion parity 

All our arguments are based on a mean-field descrip- 
tion of the proximity-induced superconductivity. Within 
this description, the number of electrons is not fixed, 
since the Hamiltonian contains pairing terms. In particu- 
lar, this means that the ground state has no well defined 
number of electrons. However, in usual superconductors 
without ground state degeneracy, the ground state |f2) 
has a well defined fermion parity. Changing the fermion 
parity requires to occupy bogoliubons, i.e. |0) — > b\ |0), 
and this always requires a finite amount of energy in a 
usual superconductor. In a TSC with ground state de- 
generacy, however, there are bogoliubons &o,m with zero 
energy. Of course, these zero energy bogoliubons also 
change the fermion parity. Thus, one half of the degen- 
erate ground states in a TSC has even parity, while the 
other half has odd parity. 

Each zero energy bogoliubon consists of two MBS, the 
wave functions of which are spatially localized. This 
means, at least on the mean-field level, that changing 
the global parity by changing the occupation of one zero 
energy bogoliubon & ,m does not affect the parts of the 
TSC that are far away from the MBS which form 60, m 
itself. The parity change affects the TSC only locally. 

In general the mean-field description of a TSC contains 
M + K bogoliubons. M of them have zero energy and 
are annihilated by the fcrmionic operators 60, m- Each 
&o, m is composed of two MBS. The remaining K ^> M 
states are finite-energy Bogoliubov excitations bk with 
energies > A. For each 60, m we define a spatial re- 
gion R rn in which the wave function of &o,m is local- 
ized. We further assume that the regions are pairwise 
non-overlapping may be specified by the local parity in 
the regions R m . As long as the TSC is in one if its 
ground states, and the positions of the regions R m do 
not change dramatically (no braiding), the local parity 
in R m cannot be changed. This is because a change of 
parity would require one bogoliubon to travel away from 
R m , and that in turn would involve extended states with 
energies > A, which are not available to the system. As 
a consequence, the topological superconductor remains 
in exactly the ground state it has been prepared in. 

B. Decoherence processes 

There are several different processes related to charge 
fluctuations on a gate that may lead to decoherence. For 
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instance, (a) a gate that is far away from any MBS may 
excite from the superconducting ground state a pair of 
high-energy bogoliubons, which travel around and acci- 
dentally hit a MBS. This is effectively similar to the phe- 
nomenon of quasiparticle poisoning liSriS Another possi- 
ble process (b) is the excitation of one single high-energy 
bogoliubon in one region R m and the departure of this 
single bogoliubon from R m . Due to parity conservation, 
this process (b) must also change the occupation of the 
zero-energy bogoliubon and thereby destroy the quantum 
information. 

While processes such as (a) are expected to depend 
on global geometric properties of the setup (how likely 
is it that an average bulk excitation finds a MBS?), pro- 
cess (b) only depends on the local conditions close to the 
MBS in question. Therefore, we focus on this process 
in this work. More explicitly, we focus on the excitation 
subprocess in region R m governed by a perturbation 
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H.c. 
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where SQ is a bosonic field that describes the charge fluc- 
tuation on the gate J£ Here we see that an excitation of a 
single high-energy bogoliubon b\. (describing an extended 
bulk state of the TSC) inevitably changes the state of the 
zero energy bogoliubon in R m . We furthermore assume 
that, once it got excited, the finite-energy excitation im- 
mediately leaves the region R m and therewith changes 
the local parity in R m . 



III. THE MODEL 

In this section we motivate the model on the basis of 
which we study the decoherence of quantum information 
via excitations of bulk bogoliubons with finite energies 
(> A). It consists of three parts. One describes the topo- 
logical superconductor -ffrsc , one describes the electrons 
on the nearby gate Hq, and one describes the coupling 
between the TSC and the gate H c . 

The TSC is modeled by an effective p-wave supercon- 
ductor in tight-binding representation^ 
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(2) 

where c„ annihilates an electron at site n, A is the super- 
conducting gap, to is the hopping amplitude, and [i n is 
the electro-chemical potential at site n. We will assume 
later that to is larger than any other energy scale in the 
TSC. In this limit, our results are independent of t,Q. In 
the proposal by Alicea et al.^ the local electric poten- 
tial on the wire is tuned via the mean charge on a nearby 
gate. However, the charge on such a metallic gate may 
fluctuate and these fluctuations may induce quasiparticle 
excitations in the TSC. Therefore, we include the gate in 
the modeling by assuming non-interacting electrons in D 



dimensions with an effective mass m* 

k 



(3) 



The operators c k annihilate an electron with momentum 
k in the gate. The capacitive coupling between the gate 
and the wire is then given by 
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with e the electron charge, e the dielectric constant of 
the insulator between the wire and the gate, eo the vac- 
uum permittivity, x„ the spatial position of the nth site 
in the TSC, r is integrated over the gate volume, and 
tp(r) = L~ D / 2 J^k e * k rc k the field operators of the gate 
electrons with L — > oo being the linear size of the gate. 
In order to make the gate-wire coupling more tractable, 
we approximate H c by assuming |r — x„| ~ d, with d the 
typical gate- wire distance, by restricting the r integral to 
a /^-dimensional Gaussian profile of width I ~ d, and by 
restricting the n summation by a profile function F„ . We 
have 




Y,F n cic n = H c . (5) 



In this work we will typically assume F n to be unity at 
the left end of the TSC (small n) where the wave function 
of the left Majorana state is non-zero, and zero on the 
right end of the wire (large n). It will turn out to be 
irrelevant how exactly F n approaches zero for large n. 
Furthermore, we may separate the mean charge (Q) from 
the charge fluctuation SQ = Q — (Q) and absorb (Q) in 
the chemical potential in -£/tsCi i- e -j we drop this term 
in H c and have 



with A = e 2 / 4ireeod. 
In the full model 
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H c excites the TSC quasiparticlcs via charge fluctuations 
5Q, the dynamics of which are governed by Hq. In or- 
der to quantify this quasiparticle excitation, it is useful 
to transform from the electron operators c„ on the TSC 
to the corresponding bogoliubon b n . For simplicity we 
restrict the discussion to the case of [i = and an odd 
number of TSC sites N in the main text and show in Ap- 
pendix [A] that these assumptions are not crucial for our 
results. After the Bogoliubov transformation one may 
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HtSC = E e k b k b k 



(8) 
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with the quasiparticlc excitation energies 

e fe = J t 2 cos 2 (fc/2) + A 2 sin 2 (fc/2) . (9) 

Note, however, that there is an additional zero-energy 
quasiparticle bo that does not appear in -Htsc- = 
lL+ilR is composed of the left/right end MBS Jl/r- Ob- 
viously, the ground state of Htsc is two-fold degenerate, 
with b bo = 0, 1. Furthermore, k = 2itm/{N + 1), m = 
l,2,..,(7V-l)/2. 

The gate- wire coupling H C7 expressed in terms of the 
Bogoliubov quasiparticles, is of the form 



\sqJ2 
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-(bkYkk' bk' 



H.c.) 
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where the sum over k, k' includes the finite-energy as well 
as the zero-energy bogoliubons. The terms containing bo 
and bk with k > describe the excitation of a local- 
ized zero-energy bogoliubon to a delocalized high-energy 
bulk bogoliubon, while terms containing only bulk bo- 
goliubons bk for k > describe the scattering between 
bulk bogoliubons or the excitation of a pair of bulk bo- 
goliubons. Since we are only interested in the processes 
that involve the Majorana state, we keep only those terms 
involving bo operators. The corresponding coupling am- 
plitudes read 
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y/N + 1 1 + S 2 - 25 cos(fc) 
with 5 = (A - t ) /(A + t Q ) and 
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t cos 2 (fc/2) + Asin 2 (fc/2) 
^t 2 cos 2 (fc/2) + A 2 sin 2 (fc/2) t 



(11) 



(12) 



Note that the equality of the modulus of the X and Y 
matrix elements for the scattering of a Majorana zero 
mode to a bulk mode or a pair creation/annihilation, 
respectively, is not special to the model or the parame- 
ters we use here, but is deeply rooted in the topological 
protection of the ground state degeneracy. It is a direct 
consequence of the fact that zero modes based on MBSs 
have equal contributions of electrons and holes in their 
wave function. 

In order to obtain a manageable expression for the 
fluctuation-induced coupling of Majorana modes and 
bulk bogoliubons [Eq. (|TT)) ] we have assumed fi = 0. As 
shown in Appendix [A~l this assumption does not affect 
our final results in an essential way. 



IV. DECAY TIME 

In this section, we estimate the time scale on which 
the system gets excited (ground state decay time) due to 



the charge fluctuations 5Q on the gate. Our final result 
has the form of Fermi's golden rule. However, due to 
the operator nature of the fluctuation 5Q, care must be 
taken in order to treat thermal and quantum fluctuations 
correctly. We therefore describe the derivation in detail. 

We assume that at time t = the topological super- 
conductor, described by Hamiltonian ((HJ), is in one of 
its ground states \ip(t = 0)) = |0). If the global par- 
ity of the TSC is even (odd), then (0|&J& |0) = (1). 
In the odd parity sector we are interested in processes 
|0) = &o|fi) — > |fc) =b\.\il), in which a zero-energy bogoli- 
ubon bo is scattered to a finite-energy bogoliubon bk- |fi) 
is the bogoliubon vacuum. The corresponding process 
in the even parity sector is |0) = |fi) — > \k) = 
where a pair of bogoliubons is created. The Hamiltonian 
restricted to these processes reads 

h = £ fe + Xok tl fc > (°l + H - c -l + Hg ' 

k fe^O 

(13) 

for both, the even and the odd parity sector. Dropping all 
terms involving matrix elements Xk,k' with both fc, fc' ^ 
corresponds to neglecting terms of higher order in per- 
turbation theory in the final Fermi's golden rule result. 

We now calculate the ground state decay rate of Hamil- 
tonian (|13p by a method based on time-dependent per- 
turbation theory. The same result can also be derived 
with non-equilibrium Green's functions on the Keldysh 
contour, but this method is significantly more involved so 
that we have chosen to describe the simpler perturbation 
theory method here. The Ansatz for the time-dependent 
state of the TSC is 



i^)) = io) + E A "E4 n) wifc) 



(14) 
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However, one must keep in mind that in this formula- 
tion the coefficient functions ai (t) are operators in the 
Hilbcrt space of the electrons on the gate, and these op- 
erators must be averaged (Og) at some point. Thus, 
the usual expression for Fermi's golden rule d|afc(i)| 2 /dt 
is ill defined. It is easily seen that an expression |afe(t)| 2 
does not respect the proper time order of the operators 
SQ [see Eqs. (fTTf and (fT5)) below]. The proper way of 
deriving Fermi's golden rule in the present context starts 
from the probability 



p k (t) = (m\k) (km)) 



(15) 



that the fcth bogoliubon is occupied. And since Pk(t) 
is itself an operator for the gate degrees of freedom, it 
finally must be averaged over the thermal occupation of 
the electrons in the gate 



Pk(t) = (Pk{t) 



Tr 



r£j. e -P(HG-e F ) 



(16) 



where j3 is the inverse electronic temperature in the gate 
and 6f is the Fermi level in the gate. The equations of 



5 



motion for the coefficients (t) are 

4 0) w = o 



(17) 



where SQ(t) = e lHat 5Qe~ lHat are the Heisenberg gate- 
charge fluctuation operators. Using the integrals of Eqs. 
(ITTl) . the leading (second) order term in Pk(t) is 



^\X 0k \ 2 dt 1 dt 2 (SQ(t 1 )SQ(t 2 )) G e- 



iefc(*l— *2) 



(18) 



for k ^ 0. In this expression, the order of the operators 
8Q(ti) and SQfa) has been properly accounted for. The 
required autocorrelation function of the gate charge 



C(t) = (SQ(t)SQ(0)) i 



(19) 



is calculated in Appendix [B] It depends on the spatial 
dimensionality D of the gate and on its linear size I. We 
find 



C d (uj) = B D cxp(-u 2 /8e F e cut ) 



[1 - exp(-0w)] 



with 



£>i — l/8ei?e cut 



B 3 = l/8e c 2 ut . 



(20) 



(21) 



Here cf is the Fermi energy in the gate, and e cu t = 
H 2 /2m*l 2 ~ 38 meV *(nm/!) 2 is an energy scale which 
must be much smaller than the Fermi energy for Eq. (|20|) 
to be valid. e cut plays the role of an ultraviolet cutoff be- 
low the band edge cutoff. It is introduced by the size I 
of the gate region which is coupled to the TSC. In the 
following discussion, the cutoff exp(— w 2 /8eF£cut) takes 
place at much higher energies than we are interested in, 
and thus it serves only as a regularization. The temper- 
ature cutoff in Eq. (|20p . however, turns out to be much 
more important. 

The spectral function Cd(uj) of the gate charge fluc- 
tuations is not symmetric in uj. For zero temperature 
one finds that C(— |w|) = 0, i.e., the quantum fluctua- 
tions are characterized by the positive frequency part of 
the spectral function. The spectral weight for negative 
frequencies is non-zero only at finite temperatures, and 
exponentially small in \oj\/T for u> < 0. 

From Eq. (jTSJ one may derive Fermi's golden rule in 
the usual way^i and one obtains for the ground state 
decay rate 



H/t = 2TrJ2\Xo m \ 2 \ 2 C D (-e m ) 

m 

= 27rB D A 2 ^|X 0m | 2 exp 



£Fe C ut/ e 



(22) 



Only the spectral weight Cd(oj) of the gate charge fluctu- 
ations with negative frequencies, i.e., the thermal charge 
fluctuations, enters the decay rate. Quantum fluctua- 
tions, which are present even at zero temperature, only 
give rise to spectral weight at positive frequencies, which 
is irrelevant for r. 

In the continuum limit (N — > oo), this rate may be 
calculated approximately in the limit to 3> A and T> A 
(see Appendix |C|) 



h 

T 



-A 2 B D Tcxp(-A/T) 



2 10 | 7 vi — (77) rcxp(-A/r) forD = l 

£ 2 £ F eV ] m \dJ *'\ I ) 

/ , \ 2 „ ,4 (23) 
281.6 ( m! ^ l[nm] nn„„„f \ /rp\ £ jj _ 3 



This is our main result. Apparently, the rate becomes 
independent of to in this limit, which is perfectly reason- 
able since to is only a model-specific detail. The physi- 
cally important parameters are the size of the SC gap A 
and the gate temperature T. 
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FIG. 1. (Color online) Decay times log 10 (r/s) as a function 
of gate temperature in Kelvin. The upper blue curve corre- 
sponds to the case A = 0.5 meV ~ 5 K and the lower green 
curve to A = 0.1 meV ~ 1 K. The dots have been calculated 
numerically from Eq. (|22[1 and the lines are the approximate 
form (|23[) . For the numerical calculation, to = 10 meV has 
been chosen. All other parameters have been set such that 
the pre-factor equals unity, TvX 2 Bn/2 = 1. 



From Fig. Q] one can see that Eq. (|23[) is a very good 
approximation of Eq. (|22[) . even for temperatures much 
lower than A. Only for A/T ~ 0.02 there are significant 
deviations from the approximate form and the pre-factor 
acquires a different power law in A/T (see Appendix |C|) . 
However, in this regime the exponential suppression is 
already so strong that the power-law pre-factor is irrele- 
vant. Thus Eq. (|23[) can be used in all relevant parameter 
regimes. 
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V. DISCUSSION 

The most important parameter entering the ground 
state decay rate r -1 in Eq. (|23]) is the ratio of the TSC 
gap A and the temperature T. For T <C A the rate t -1 
is exponentially suppressed by a factor exp(— A/T), as 
expected for a system in thermal equilibrium. However, 
the pre-factor gives hints about the optimal gate design. 
First of all one can see that there are huge differences 
between a ID and a 3D gate. This is due to the different 
pre-factors Bp in the correlation function of the gate 
charge fluctuations. For realistic parameters one finds 
B3/B1 ~ 10 3 . For a ID gate, the pre-factor is of order 
one for / = d (linear size I of the active gate region is 
equal to the gate- wire distance), e = 1, cf = lcV, and 
to* = O.Itoo. For a 3D gate, the gate size I does not scale 
with the gate- wire distance d. With the same parameters 
and I = d = 50 nm the pre-factor is ~ 7000. From 
this one may draw the conclusion that the optimal gate 
has (effectively) low dimensions, a small active region, 
and is not too close to the TSC. Similarly, the inverse 
dependence on the dielectric constant e and on the gate 
material parameters to* and €f gives further hints for 
the optimal gate design. 

We are now in a position to give reasonable limits for 
the two important parameters A and T. For this we 
calculate the minimum gap A m j n ~ T\og(Tir\ 2 B]jT/2h), 
required for a topological protection that should last for 
at least one microsecond, i.e. r > 1 (is. Of course, A m j n 
depends on the pre-factor ttX 2 Bjj/2, which may be tuned 
by several orders of magnitude. We assume here a worst 
case scenario ttX 2 Bd/2 = 10 4 for a 3D gate with I = d = 
50 nm, and a best case scenario tt\ 2 Bd/2 = 0.01 for an 
optimized ID gate (e.g., dielectric constant e = 5 and 
l/d = 1/2). Figure [2] shows A m i n for gate temperatures 
smaller than 1 Kelvin and for the best and worst case 
scenario. 

Approximately, a TSC gap A > 20 T is sufficient to 
obtain a TSC, the ground state of which is stable for 
at least one microsecond. By optimizations in the gate 
design this ratio may be reduced by a factor ~ 4. 

VI. CONCLUSION 

We have studied the loss of quantum information 
stored in a Majorana qubit due to thermal charge fluctu- 
ations on a nearby gate, which is capacitively coupled to 
the topological superconductor (TSC) hosting the Ma- 
jorana fermions. The fluctuating gate excites the TSC 
above its ground state. Since the topological protection 
is only active as long as the TSC is in its ground state, the 
stored information gets lost as soon as a finite-energy ex- 
citation is created. Starting from a tight-binding model 
of a TSC,— we derived a microscopic description of the 
gate-TSC coupling via the Coulomb interaction. Finally, 
we obtained an expression for the qubit lifetime which 
only depends on the gate temperature T, the supercon- 
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FIG. 2. (Color online) Minimal gap of the TSC A min in Kelvin 
as a function of the gate temperature T. The upper green 
(solid) line corresponds to the worst-case scenario in which the 
rate pre-factor ■k\ 2 Bb/2 — 10 4 , while the lower blue (dashed) 
line refers to the best-case scenario for which n\ 2 B D /2 ~ 0.01 
(see the main text). 



ducting gap A, and some parameters describing the gate 
geometry and material parameters. Our result does prac- 
tically not depend on the hopping amplitude and the 
chemical potential of the TSC model. 

It was found that A/T ~ 20 is required for a qubit 
lifetime of at least one microsecond, if a non-optimized 
gate design is chosen. By using optimized gate geometries 
and dielectrics, however, this ratio may be reduced by 
factors of 4 and more. 

This work shows that not only the temperature of 
the topological superconductor, in which quantum in- 
formation is stored, is relevant for the coherence time. 
Also the temperature of distant parts in the experimental 
setup, coupled to the superconductor via the long-ranged 
Coulomb interaction, is indeed important and must be 
sufficiently small compared to the gap of the supercon- 
ductor. On top of that, we mention again that any exper- 
imental realization of a topological superconductor still 
suffers of the seemingly ubiquitous problem of quasipar- 
ticle poisoning. 

Finally we note that thermal fluctuations are not the 
only relevant source of gate noise. For a functional quan- 
tum computer based on Majorana qubits it is also im- 
portant to investigate the non-equilibrium noise due to 
the changes in the gate voltages, needed for the actual 
braiding operations. This issue may be studied on the 
basis of the model we have derived, but is beyond the 
scope of this work. 
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The corresponding energies are given in Eq. ([9]) and are 
equal for both subchains. 

Knowing the Bogoliubov transformation matrices, 
which can be read off Eqs. (|A3|) - (|A5|) . one may rewrite 
the gate- wire coupling H c in terms of bi [see Eq. (fTU|) ]. 
Straightforward algebra shows that the matrix elements 
Xqi of H c are given by 



Appendix A: Gate-wire coupling for /J. ^ 

In the main text we have restricted the discussion to 
/i = in order to obtain manageable expressions for the 
gate-wire coupling. In this section we show that relaxing 
this assumption does not change the results dramatically. 

The Bogoliubov transformation from the electron basis 
c„ to the bogoliubon basis bk in the topological supercon- 
ductor reads 



(Al) 



I is not necessarily related to a momentum but should 
be thought of as a general label of the Bogoliubov eigen- 
states with energy e;. n is a real space index, labeling 
the sites of the TSC wire. For fi = and TV odd the 
Bogoliubov transformation becomes particularly simple 
because -f/rsc can be decomposed into two uncoupled 
chains with alternating hoppings when written formally 
in terms of Majorana operators 



-ffrsc = i 



(A + t )"fBnJAn 



+ 1 



+ (A - to)~fAnlBn+l - Z^lAnlBr 



(A2) 



where c„ = ja,u + i"fB,n- Moreover, because of N odd, 
the two end Majoranas are located on different subchains 
and can therefore be treated separately even for finite N. 
For the zero-energy bogoliubon one finds 

(N-l)/2 



K = ^ ^ n (C2n+l+ci n _ 1 + cJ v _ 2n -CjV-2n), 

n=0 

(A3) 

with S = (A — 1 ) /(A + t ) and for the bulk bogoliubons, 
where the index is I = (k,s) with s = 1,2 for the two 
subchains, we have 



'k,l 



(N+l)/2 

E 



n=l 



sin(fcn)(c 2 „ - c\ n ) 



sin(fcn + al)(c 2n -i + 4„-i) 



(A4) 



{N+l)/2 

E 



n=l 



sm(kn)(c 2n + 4n) 



+ sin(fcn + al)(c 2n -i - 4n-i) 



(A5) 



AV = ai, n F n a^ n - fr, n F n ft, in (A6) 
Y lM = Pt n F n at, !n - al n F n ft, tn . (A7) 

If one of the indices I or V corresponds to the zero mode, 
and this is the only case we are interested in, then one 
can easily show that |A"o./| = \Yq,i\ = l^i.ol = l^i.ol an d 
it is thus sufficient to consider Xq.i- Furthermore, since 
OiQ <n and /?o, n are exponentially small for n away from 
one of the ends, we set F n = 1 for n < N/2 and zero 
otherwise. This corresponds to a gate which affects the 
TSC over a region which is at least as large as the spatial 
size of the MBS wave function. The latter is typically 
on the order of 100 nm. Since the wave function of the 
MBS (ao, n and Po,n) is exponentially small for n ~ N /2 
where F n deviates from unity, the particular form of the 
profile function F n near N/2 is not important and F n can 
be chosen to be technically most convenient. Evaluating 
Eq. (|A6|) under those assumptions (i.e., /i = 0; N odd; 
F n = 1 for n < N/2 and otherwise) results in the simple 
analytical expression Eq. (fTTT) . Note that the Xqj is only 
non-zero if I = (k, 2), i.e., for the second subchain. 

In the remainder of this section we calculate the Bogoli- 
ubov transformation ai t „ and numerically from the 
eigenvectors of the TSC wire Hamiltonian ifrsc m Eq. 
(P| with an arbitrary number N of sites and for general 
/j,. For (i 7^ the subchains arc coupled. From ai. n and 
Pi,n we calculate Xq^ numerically and use it to evaluate 
the Fermi's golden rule expression for the scattering rate 
1/t [Eq. i]22p]. The main effect of a nonzero chemical po- 
tential on the scattering rate is a trivial renormalization 
of the superconducting gap 



a*(aO = aJi- 



(A8) 



Thus, in order to resolve the nontrivial correction, we cal- 
culate the scattering rate for /i-dependent temperatures 
T = A*([x)/M, with M = 1,10,100, and compare the 
result to the scattering rate at fi = 0. Figure [3] shows 
that t depends only weakly on the chemical potential. 
Typically, the scattering rate is renormalized by factors 
between \ and | for > 0. Note also that the chemical 
potential is varied over a relatively large range ~ to in 
Fig.H 

Furthermore, we have calculated Fig. [3] for N = 500 
and iV = 501. The difference is not visible. Thus, the 
simplifying assumptions N odd and /i — are well justi- 
fied. 
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FIG. 3. (Color online) The ratio r(fi = 0)/r(/i) of the 
scattering rates at different chemical potentials /i. We 
used the parameters A = 0.1, 0.2, 0.6, to — 1, T = 
A*(aO,A*(aO/10, A*(ju)/100, € F = 4 and e cut = 0.004. The 
number of sites in the TSC wire in this calculation is N = 501. 
The result of a calculation with TV = 500 cannot be distin- 
guished from this figure. The curves with larger deviations 
from 1 correspond to larger A and smaller temperatures. 



Appendix B: Gate correlation functions 



In this appendix we outline the calculation of the cor- 
relation function of the gate charge fluctuations C(t) = 
(SQ(t)SQ(0)) G , where 



Tr [e 



-p{H G -e F ) . , 



(Bl) 



In frequency space and D dimensions we have 

C D ( U ) = (J^\ j d^kd^qe-^/ 2 ^ 

S(u - E k + £ k+q )n k+q (l-n k ), (B2) 
I 



with the Fermi factors 

n k = (expC3(£ k - e F )) + l)" 1 , 
and the electronic energy in the gates 

a " 2 



2m* 



Furthermore, the cutoff energy scale 

1 



Ccut 



2m* I 2 



(B3) 



(B4) 



(B5) 



has been introduced. We now calculate the integral (|B2 
for one and three spatial dimensions. 



1. D = 1 

We first eliminate the k integral via the Dirac delta 
function and transform the q integration to an energy 
integration E = E q . This gives 



Ci(w) = 



1 



AE 



-B/2e cu 



87re cut Jo E 
with the smeared plateau function 

V{w,x,T) 



0J 

AE 



E 



(B6) 



,-{x+w/2)/T + 1 yl, e (x-w/2)/T + ^-1 

(B7) 



describing a plateau with height 1 and width w centered 
around x = 0. The plateau edges are smeared over Ax ~ 
T, 

The second argument of the plateau function in Eq. 
(|B6[) is close to zero for E\ = Aep and for E% = uj 2 /4cf, 
where we have assumed that u < 6y. Because of the 
exponential cutoff e cut <^ e F we may neglect the first 
energy region. We do this by pulling the exponential 
function out of the integral at E = Ei and at the same 
time regularizing the integral such that the integrand is 
exponentially small for energies near E±. The simplest 
way to do this is to change the sign of the E/A term in 
Eq. (|B6[) . This does not affect the integral for energies 
around E^ but changes the sign of E\. Thus, we have 



1 



87re C ut 
1 

87re cu t 
1 

87re cut 
1 



exp 



cxp 



exp 



8e F e 



cut / Jo 



dE 
~E' 



^ E 

u), e F , T 

' 4E 4 



uj 2 \ r 2dy e"/ 2T 



2dy 

^F^nt) J-oo E + uj 2 /E cosh(w/2T) + cosh((t/ - e F )/T) 

,2 \ 2 r°° p lj / 2T 



8ire F e c 



■ exp 



•(-F^cutJ 4e^ 
oj' 



-T- / dy 



cosh(o;/2T) + cosh((y - e F )/T) 



,2 \ 



8eFe cu t J I - exp(-w/T) ' 



(B8) 
(B9) 
(B10) 
(Bll) 
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where we have introduced y = to 2 /AE — E/A. Further- 
more we have used that the integrand is only large for 
E ~ E-2- All approximations become exact in the limit 
ep — > oo. For typical ep ~ eV the relative errors of the 
approximations used are on the order of 10 -5 . 



D 



For a gate with three spatial dimensions we first elim- 
inate the integral over the angle between the vectors k 
and q via the Dirac delta function and obtain 



C 3 (w) 



8tt 



dkkn(E k -w)(l -n{E k )) 



'■/+ 



dq qe 



-lV/2 



(B12) 



with fc_ = Re V 2mw and g± = |fc ± \Jk 2 + 2muj\. With 
y = k /2m* — lu/2 one may write 

C 3 (w) = 



logout 



50 , cxp(- ? //e cut ) sinh(^y 2 - ^ 2 /4/e cut ) 

, /2 , y cosh(w/2T) + cosh((y — ep)/T) ' { 6> 



For ep T,uj,e cu t this integral may be approximate as 
described above in the D = 1 case and one obtains 



c 3 H 



1 



■ cxp 



8e F e CVL tJ 1 - exp(-w/T) 

(B14) 



Appendix C: Approximative evaluation of the 
excitation rate 

In this appendix we outline the approximation of the 
Fermi's golden rule expression of the excitation rate given 
in Eq. (|22[) . In the limit ./V — > 00 the sum over m becomes 
a k integral and r _1 is proportional to 



(l-<5 2 ) / dfc 
'0 



(sin(fc + a k ) - Ssma k ) 2 e k 



(1 + 5 2 - 25 cos A:) 2 



(CI) 

After transforming the integral to z = to(-K — k)/A, the 
limit to/A—> co may be performed for /. We define 



I* = lim / 

tn— ¥00 



= 32A 



{4 + z 2 )^(exp(£VI+^)-l) 



In realistic situations, to may be assumed to be at least 
two orders of magnitude larger than the induced su- 
perconducting gap A. In this case, the relative error 
(/ — I*) /I is on the order of 10~ 4 . Thus, using /* is an 




-0.4 



0.00 



0.05 



0.10 0.15 

T/A 



0.20 



FIG. 4. (Color online) Relative errors (1 — / a pprox/-0 for var- 
ious approximate expressions / ap prox of the rate integral / as 
a function of T/A. Part (b) is a zoom of part (a) for small 
values of T/A. 



excellent approximation. In this way, the hopping pa- 
rameter to of the TSC model, which has not much direct 
physical meaning, is removed from the rate. 

The expression (|C2[) is dominated by the exponential 
factor exp(— A/T). We thus continue by evaluating the 
A/T -> limit of re A ' T 

/•OO 2 rp 

lim Pe A / T = 6AT / dz - ... = ^f, (C3) 
a/t^o J (4 + z 2 ) 3 2 v ' 



(C2) go that for A«T 



/•~^exp(-A/T). 



(C4) 



Note that A/T = is a singular point of I* so that the 
expansion in A/T can be performed up to the fourth 
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order. However, this expansion has zero convergence ra- 
dius. It becomes better for A/T < 10 up to second order 
and then becomes worse again. In second order one finds 



I* 



A 

T 



5tt 
12 



(C5) 

For A/T > 10 this approximation is much worse than 
the simple limit (|C4|) . 



Next, we evaluate I* in the limit of large A/T. In this 
case the exponential function is much larger than unity 
so that the 1 may be neglected in Eq. (|C2[) . Furthermore 
one may expand the square root in the argument. One 



finds 

r°° ? 2j a 
r= i 32Ajf I 4 T ^ F ,exp(--(l + z 2 /8)) 

= 2^(3/2; ;A)Ae- A / r , (C6) 

with Tricomi's confluent hypergeometric function 
C/(a;6;2). From Fig. 2] one can see that this approxi- 
mation becomes worse than (|C4[) for T/A > 1/2. The 
dashed line in part (b) of Fig. [4] shows an expansion of 



[7(3/2; b; z) up to order (T/A) 7 / 2 



i.e. 



'5 



'I 



Ae- A / T . 



(C7) 

This expansion is only valid for extremely low temper- 
atures where the exponential suppression dominates the 
rate. 

Thus we may conclude that, although (|C6[) is a very 
good approximation for T/A < 0.12 and JCS) for T/A > 
0.12, the simplest form (|CJ4|) is sufficiently accurate in the 
relevant range down to T/A ~ 0.02. Below this temper- 
ature the exponential suppression of the rate is so strong 
that the pre-factor is irrelevant. 
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